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Abstract: This paper is concerned with a networked estimation problem in which sensor 
data are transmitted over the network. In the event-based sampling scheme known as level- 
crossing or send-on-delta (SOD), sensor data are transmitted to the estimator node if the 
difference between the current sensor value and the last transmitted one is greater than a 
given threshold. Event-based sampling has been shown to be more efficient than the time- 
triggered one in some situations, especially in network bandwidth improvement. However, 
it cannot detect packet dropout situations because data transmission and reception do not 
use a periodical time-stamp mechanism as found in time-triggered sampling systems. 
Motivated by this issue, we propose a modified event-based sampling scheme called 
modified SOD in which sensor data are sent when either the change of sensor output 
exceeds a given threshold or the time elapses more than a given interval. Through 
simulation results, we show that the proposed modified SOD sampling significantly 
improves estimation performance when packet dropouts happen. 

Keywords: Networked estimation; event-based sampling; send-on-delta; packet dropout. 



1. Introduction 

Recent works have discussed event-driven alternatives to traditional time-triggered sampling 
schemes. It has been shown to be more efficient than time-triggered one in some situations, especially 
in network bandwidth improvement. In [1-7], event-based sampling scheme was applied by adjusting 
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the threshold value at each sensor node, data transmission rate is reduced so that the network can be 
used for other traffic. 

However, analysis and simulation in the the works on event-driven sampling scheme were 
performed under ideal communication network conditions: no delays or packet dropouts are assumed, 
but in realistic applications, network induced delays and packet losses do happen. 

The issues of network delays and packet dropouts in time-triggered systems have been addressed 
and solved by researchers in [8-14]. In [8] the stability of the Kalman filter in relation to the data 
arrival rate is investigated. It is shown that there exists a critical data arrival rate for an unstable system 
so that the mean filtering error covariance will be bounded for any initial condition. In a very recent 
study [13], the optimal H 2 filtering problems associated respectively with possible delay of one 
sampling period, uncertain observations and multiple packet dropouts are studied under a unified 
framework. The Eb-norm of systems with stochastic parameters is defined and computed via a 
Lyapunov equation and a steady-state filter is designed via an LMI approach. In [14], the authors adopt 
a model similar to that of [13] for multiple packet dropouts to investigate finite-horizon optimal linear 
filtering, prediction and smoothing problems. 

In conventional event-based sampling systems, also called send-on-delta (SOD) sampling [5-7], the 
issues of network delay and packet loss are difficult to solve because data transmission and reception 
do not use a periodical time-stamp mechanism as in the time-triggered sampling systems. Motivated by 
those issues, in this paper, we introduce a modified SOD sampling scheme in which the event-driven 
sampling is combined with a time-triggered sampling scheme to detect packet dropouts. Then, a 
networked estimator based on a Kalman filter is formulated to estimate states of the system 
periodically even when the sensor nodes do not transmit data. The proposed SOD sampling scheme 
has properties inherited from the conventional SOD sampling: so the benefits from event-driven 
sampling are still hold. Through theoretical analysis and simulation results, we show that the proposed 
SOD sampling scheme gives better estimation performance than the conventional SOD one when 
packet loss happens. 

2. Modified SOD Sampling Scheme 

Consider a networked control system described by the linear continuous-time model: 



where x(t) € R n is the state of the plant, u is the deterministic input signal, y(t) G R p is the 
measurement output which is sent to the estimator node by the sensor nodes. w(t) is the process noise 
with covariance Q , and v(t) is the measurement noise with covariance R . We assume that w(t) and 
v(t) are uncorrected, zero mean white Gaussian random processes. 

The modified SOD sampling scheme illustrated in Figure lb is stated as follows: 
Let yi astj i (1 < i < p ) be the last transmitted value of the z-th sensor output at instant t last ^ . A new 
sensor value will be sent to the estimator node if one of two following conditions is satisfied: 



x(t) = Ax(t) + Bu(t) + w(t) 

y(t) = Cx(t) + v(t) 



(i) 



Vi(t) - Vlast,i \ > s t 



(2a) 



Sensors 2009, 9 



3080 



t - thut,i > 6t,i ( 2b ) 
where 6 y , n 6 t i are the given magnitude, time threshold values respectively at the z-th sensor node. 

Figure 1. Principle of SOD and modified SOD sampling schemes. 
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Using the modified SOD sampling scheme above we will obtain some benefits. Firstly, the 
estimator can detect signal oscillations or steady-state error if the difference of output value remains 
within the threshold range during a long time. Secondly, the estimator can detect multiple packet 
dropouts if it does not receive sensor data within the interval (0,8 t i). Thirdly, theoretical analysis for 
SOD sampling is still applied for the modified SOD sampling. 

However, this scheme has one disadvantage that sensor data transmission rate will be increased due 
to condition (2b). If 6 ti is small, the estimator detects packet dropouts fast but data transmission rate is 
increased. If 6 ti is large, transmission rate is small but the estimator detects packet dropouts slowly. 
Therefore, an optimal 6 t j value is necessary to compromise these constraints. 



2.1. Multiple packet dropouts detection 



The estimator node detects packet dropouts of z'-th sensor data by checking the instant z'-th sensor 
data arrive. If there is no z-th sensor data arriving, the estimator node for the time t — ti ast l > 6 tt , then 

the estimator node knows that one-packet dropout happened at the z'-th sensor node. Similarly, if there 
is no z'-th sensor data arriving for t — ti ast . t > 28 t , l , then two-consecutive-packet dropout happened. 

We state the general case for multiple packet dropouts as follows: 

If the estimator node does not receive i-th sensor data for time (t — ti astii ) > diS ti (d { = 1,2,3,...) 
then the estimator knows that there have been at least d consecutive packet dropouts at the i-th sensor 
node since the time receiving yi ast) i ■ 

Note that the estimator just detects "at least" d { consecutive packet dropouts, not precise d { 

consecutive packet dropouts because there exists a delay interval in detecting packet dropouts. As 
illustrated in Figure 2, although packet loss happens within the time range (ti ast)i: ti ast)i + St,i)> the 
estimator only detects it at a time (ti ast , + 6 t , t ) . Thus, if there is more than one packet dropout within 
the time range (ti ast) i : t last) i +6 ti ), the estimator also detects only one packet dropout at time 
{kast,i + 8t,i) ■ This is an inevitable flaw of the modified SOD sampling scheme. We can constraint this 
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flaw by reducing the 6 ti value, but sensor data transmission rate will be increased. Therefore, an 
optimal 6 ti value is necessary to compromise between the two constraints. 

Figure 2. Multiple packet dropout detection. 



Hi W Packet dropout At time the estimator detects 
happens packet losses 




8 t 8 t 6 t 8, S t 



3. State Estimation with Modified SOD Transmission Method 

The networked estimation problem applying modified SOD transmission method can be described 
as follows: 

1 . Measurement output y { ( 1 < i < p ) are sampled at the period T but their data are only sent to 

the estimator node when (2a) or (2b) is satisfied. 

2. For simplicity in the problem formulation, transmission delay from the sensor nodes to the 
estimator node is ignored. 

3. The estimator node estimates states of the plant regularly at the period T regardless of whether 
or not sensor data arrive. If there is no i -th sensor data received for ( t — ti ast 4 ) > d^ , the 
estimator node considers that the measurement value of the i -th sensor output y i (t) is still equal 
to yiast j but the measurement noise increases from v^t) to v ni (t) = v t (t) + Aj(^,^ ast ,) . 

Note that if = 0 then there is no packet dropout, the estimator acts like a conventional SOD 
filter [5]. To formulate a state estimation problem, the boundry of A.i(t,ti ast $) needs to be determined 
as di ^ 0 (packet dropouts happen). In the next section, we will compute the covariance of v nji (t) 
when d i ^ 0 and then a modified Kalman filter is applied for state estimation. 

3.1. Measurement noise increased due to multiple packet dropouts 

We know from (2a) that ly^t) — y[ asti \ < S y i as long as the estimator node does not receive a new 

i -th sensor data value. If one packet dropout happens, the i -th sensor output value has changed more 
than 6 y i . The estimator should know that: 

\Vt(t) ~ Vlasu \ < 8y,i +$y,i 
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For general cases, as shown in Figure 3, if there are d { consecutive packet dropouts then: 

A i(t,t h st,i) = \Vi(t) ~ Vkut,i\ < ( d i + tyy.i ■ 



(3) 



Note that (3) is also applied to the case of no packet dropout [5] by letting d t = 0 . Assuming that 
A,(£, ti astji ) has a uniform distribution with (3), variance of A t (t : ti as ti) will be: 



(4) 



£[A;(M, asM )] = 0 
E[A^t,t last!i )] = ((d i+ l)6j/3 
Var[A,(t,t las j} = E[^(t,t las j] - E 2 [A^J 
= ((d, +l)6 yA f /3 

Therefore, if there is no i -th sensor data received for t > ti ast ^ , variance of measurement noise is 
increased from R(i, i) to R(i, i) + ((di + l)6 yi f / 3 . 

Figure 3. Measurement noise increased due to multiple packet dropouts. 
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3.2. State estimation 

A modified Kalman filter for state estimation x k at step k, where there is a change in the 
measurement update part of the discrete Kalman filter algorithm [15], is given as in the Figure 4. We 
use the discretized system model sampled at period T : 

T 

A.=e AT , B d = J e Ar Bdr, 
o 

where Q d is the process noise covariance of the discretized system: 



Q d = fe Ar Qe A ' r dr, 



and y last is the vector of p last received sensor values: 



Via 



Ulastl Ulast.2 



Ulast.p 
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Figure 4. Structure of the modified Kalman filter. 
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X k+1 — Al X k + B d U k 
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In the modified Kalman filter in Figure 4, the states of the plant are estimated regularly at every 
period T, regardless of whether or not sensor data arrive. If i -th sensor data arrive then 
A-i(t,ti astji ) = 0, the modified Kalman filter acts like the conventional Kalman filter. Otherwise, if i- 
th sensor data do not arrive due to packet loss, it uses yi ast ,i as the measurement value and 
R(i,i) = R(i,i) + ((d, + l)6 y i f / 3 as measurement noise covariance for state estimation. 

As stated in [8], if the system (1) is unstable and a packet loss rate is high, the proposed filter could 
diverge. For example, if all packets are lost, d t will increase and thus R t will become infinite. Thus 
Pin Figure 4 could become infinite. 

4. Optimal ^ Computing Problem 

As mentioned in Section 3, S t ,i is a trade-off parameter between sensor data transmission rate and 
the response of packet dropouts detection. The response of packet dropout detection guarantees 
estimation performance. Because SOD sampling is more efficient than the time-triggered one in 
network bandwidth improvement, we should choose d t j such that sensor data transmission rate is 
reduced to promote ability of SOD sampling. In the next section, we will investigate the relation of d t ,i 
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with transmission rate and the effect of d tii on estimation performance. Then an optimization problem 
is formulated to find the optimal d tii value according to the given estimation performance. 

4.1. Sensor data transmission rate by condition (2b) 

The total sensor data transmission rate caused by condition (2b) in a time unit: 



i=l °t,i 



(5) 



where p is the number of sensor output 



4.2. Estimation error covariance due to packet dropouts 



Let £j (0 < £j < 1) be the packet loss rate at the z'-th sensor node, 4; = 0 corresponds to no packet 
loss. Let A T { be the average transmitting time per packet of the z'-th sensor node in the conventional 
SOD method. Note that A T { is dependent on the given 8 Vti value, but independent on d t j value. A T { 

is computed by running the simulation model in analysis. In practice, it can be computed by letting 
6 ti = oo and monitoring the number of packets in a time unit. 

The average number of packet dropouts in the conventional SOD sampling per a time unit: 

6 



1 AT 



(6) 



In the proposed SOD sampling, the average number of packet dropouts within the time interval 6 tii 
will be: 

Mi (7) 



cL 



AT 



We know from Section 4.1 that the larger number of consecutive packet dropouts is, the larger 
measurement noise covariance is. Measurement noise covariance is largest if di packets are 
consecutively lost. Following the idea in (4), if there is d t packet loss, the measurement noise 
covariance should be increased as follows: 

\2 



= %o + 



AT 



+ 1 



(8) 



4.3. Optimal 6 ti computation 



In this section, d t j value is computed. Using (8), we assume that the measurement noise covariance 
is given by: 



R = R + Diag 



[AT, + 



AT 



2 



(9) 
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The estimation performance in this case can be computed from the following discrete algebraic 
Riccati equation: 

P = A d PA' d +Q d - ^'(CPC + RY'CPA, (10) 

Note that (10) does not provide the actual estimation error covariance of the filter. The main 
purpose of (10) is to evaluate how d tii affects the estimation performance. We can see that if S t j is large, 
the estimation error covariance P increases. 

The solution of (10) is denoted hyP(8 t j ) . In the following optimization algorithm to find Si, we try 

to reduce the sensor transmission rate caused by condition (2b) subject to the given estimation 
performance constraint: 

d t Optimization Problem 

(11) 

subject to DiagP(8 t i ) < fj,P 0 

where P 0 is the upper bound error covariance with given value 6 y>i and no packet dropout (solution of 
(10) as d = Diag(0,...,0)). P 0 is also the estimation performance of the conventional SOD. \x is the 
ratio to the estimation performance of conventional SOD filter in case of no packet dropout. If \i is 
large, the 6 t i optimization problem (1 1) is done with weaker estimation performance constraints. 



5. Simulation 



To verify the proposed filter, we consider an example of the second-order system with step input 
where the output is sampled by the SOD and modified SOD sampling: 





0 






0 




x(t) = 


-1/a 


-b/a 


x(t) + 


M/a 


u(t) + w(t) 



y(t) 



1 o 



30, a = 5, b = 1 



x(t) + v(t) 
Q = 0.01, R = 0.01, T = 10ms 

where the system parameters for performance evaluation are given by M 

(underdamped system) . The simulation process is implemented for 50 seconds. 

Choose n = 5 for the optimization problem (11). The solution 6 t i, 6 t 2 of (11) along with 8 y i and 

are shown in Figures 5 and 6, respectively. We see that 6 tji is proportional to 6 y j and reversely 

proportional to £ 4 . It means that when 6 y j is large, the z'-th sensor data transmission rate is small, thus 

6 ti is also small to keep the overall transmission rate small. But if packet dropouts increase ( & is 

large), 8 ti value is lowered. As the result, the overall sensor data transmission rate is increased to 

guarantee estimation performance. 
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Figure 6. <5 t 2 of (11) along with 6 y2 and £ 2 




Table 1. Estimation error along with packet loss rate in two filters. 



Packet loss rate 

& = & 


0.05(5%) 


0.1(10%) 


0.15(15%) 


0.2(20%) 


n (SOD) 






ni 

n 2 


= 95 
= 31 








4, 


t =4.12 


8 tl =2.08 


8 tl = 1.73 


k 


i = 1.52 


ki 


4. 


2 =4.69 


6 t2 =2.31 


6 t2 =1.91 


k 


2 = 1-66 


n (modified SOD) 


ni 

n 2 


= 101 

= 36 


m = 109 

n 2 = 44 


m = 112 

n 2 = 47 


ni 

n 2 


= 115 
= 50 


e (SOD) 


ei 


= 0.0383 


ei = 0.0384 


ei = 0.0386 


ei 


= 0.0391 


e 2 


= 0.0167 


e 2 = 0.0168 


e 2 = 0.0169 


e 2 


= 0.0172 


e (modified SOD) 


ei 
e 2 


= 0.0075 
= 0.0096 


ei = 0.0064 
e 2 = 0.0089 


ei = 0.0039 
e 2 = 0.0082 


ei 
e 2 


= 0.0020 
= 0.0069 



Table 1 shows the estimation error in two filters (SOD filter and modified SOD filter) as 
k,i = = 0-5, /i = 5 and ^,^ 2 are varying 5%, 10%, 15%, 20%. Estimation error is evaluated 

by: 
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1 N 



N 



„ \2 
■Ek,i ) 



(12) 



fe=i 



where x, is the reference state, x t is the estimated state, and N = 5,000. 

In Table 1, we see that when applying the modified SOD filter, the estimation error is significantly 
improved. For instance, in the case & = £ 2 = 0.05 , the total number of sensor data transmissions in 

the modified SOD (# 137) is just slightly greater than that in conventional SOD (# 126) but the 
estimation error is reduced so much ((ei = 0.0075, Q2 = 0.0096) compared to (ei = 0.0383, Q2 = 
0.0167)). 



Figure 7. Estimation error in two filters as £ x = £ 2 
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Figure 8. Instants the sensor node transmits data due to condition (2b). 
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Figure 7 intuitively shows the estimation error in two filters as & = £ 2 = 0.05 , 8 y i = 6 y 2 = 0.5, 
6 t i =4.12, 6 t 2 = 4.69 . The boundry of e x in the modified SOD filter (SODa) is much smaller than 
that in the conventional SOD filter. Figure 8 shows the instants the sensor node transmits data to the 
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estimator node due to condition (2b). We see that the number of sensor data transmissions caused by 
condition (2b) is very small in comparison with the total number of sensor data transmissions [(ni = 7, 
n 2 = 7) compared to (m = 101, n 2 = 36)]. When the modified SOD sampling is applied, the total 
number of sensor data transmissions is slightly increased, but the estimation error is significantly 
reduced. Therefore, the modified SOD sampling significantly improves estimation performance with 
only a little increase in the data transmission rate. 

Notice that if we just consider the transmission condition (Equation 2a), estimation error of the 
proposed method is worse for systems that the output varies slowly. However, an issue of conventional 
event-based sampling is that it can not detect signal oscillations or steady-state error if the difference 
of output value remains within the threshold range (because the output varies slowly). This fact causes 
estimation error to be increased. Whereas, the proposed method uses the transmission condition 
(Equation 2b) not only to detect packet dropouts but to reduce the error in case the output changes 
slowly. 

As illustrated in Figures 7 and 8, where the estimation error of the proposed method (top-right graph 
of Figure 7) and of the conventional method (top-left graph of Figure 7) are shown according to the 
output yl (top-left graph of Figure 8). We see when yl varies slowly (time interval from 20s to 50s), 
the proposed method gives much smaller estimation error than the conventional one. 

In case the output changes fast, it is obvious that ignoring packet dropout will introduce extremely 
incorrect result because we still use the wrong old measurement noise value even when we do not 
know how much the output value changes. 

5. Conclusions 

In this paper, the state estimation problem with modified SOD transmission method over networks, 
in which an event-based sampling is combined with a time-triggered sampling to detect packet loss 
situations, has been considered. We have shown that when using the proposed modified SOD filter, 
estimation performance is significantly improved with a small increase in sensor data transmission. If 
multiple packet dropouts happen, the estimator node will detect and compensate for them with an 
amount of additive measurement noise to improve estimation performance. This method is very useful 
for networks where data transmission is unreliable due to noise. 
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